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Abstract 

Using the symmetry approach, we find a class of integrable nonlinear PDEs with 
dispersion law u>(k) = k^. All these equations turn out to be linearizable by means 
of a differential parametrization. 

1 Introduction 

We consider the problem of classification of PDEs 

qtt = Qxxx + F(q, q x , q t , q xx , q tx ). (1) 

having higher symmetries. These equations possess an unusual dispersion law u(k) = k*. 
Moreover one can prove that for such equations there are no higher conserved densities 
and therefore they are nonhamiltonian. 
The first example 

qu = q X xx + 3q x qtx + {qt - 3q 2 x )q xx (2) 

of such equations has been found in [1]. Particular solutions of (2) feature moving branch- 
points singularities. For instance, the simplest self-similar solutions of the form q = 
q(x + kt) satisfy the ODE 

y' — y(y — k) 2 + const 

where q x = y. It is obvious that they have singularities of such type. 

Nevertheless, equation (2) possesses a degenerate Lax representation which implies 
infinitely many higher symmetries and corresponding invariant solutions [1] . Three higher 
symmetries are as follows 

q tl = q xx + 2q x q t - q\ 

Ota = qxt + q x q xx + qt + q\qt - ql 

q t3 = qtt + 3q x qt ~ 2q%qt- 
The whole hierarchy of symmetries can be reproduced by a recursion operator 

TZ = -q x + 2D' l q xx + D~ x D t = A + BD t (3) 
acting on the seed symmetries q x and q t . 



In Section 2 we generalize the main concepts of the symmetry approach, such as the for- 
mal recursion operator and canonical conserved densities for the case of non-evolutionary 
equations of the form 

Qtt = F(q, qi, q 2 , . . . , q n , q t , Qti, Qa, ■ • • , Qtm)- (4) 

Such type of equations were excluded from consideration in works [2, 3, 4], where only 
evolution equations were investigated. Obviously, equation (4) can be rewritten as a 
system of two evolution equations. For example, equation (2) is equivalent to 

iH = {v x -u v) x , v t = (u- v 2 ) x , 

where v = q x . However the matrix coefficients of leading derivatives for such systems 
have the structure of a Jordan block, whereas in papers [5, 6, 7], devoted to systems of 
evolution equations, the leading matrix was supposed to be diagonalizable. 

Using the technique of canonical conserved densities, we find all equations (1) possess- 
ing infinitely many higher symmetries. It turns out (see Section 3) that all these equations 
are related to second order evolution equations 

iH = H(x,u,u x ,u xx ) (5) 

having higher symmetries. It was shown by Sergey Svinolupov [8] that any such equa- 
tion (5) can by reduced to one of the following equations 

u t = u xx + f(x)u, 

u t = u xx + 2uu x + g(x), 

u t = ("I + Xx ) ' 

\ u A / x 

(U x \ 

u t — [ + X\xu + X 2 u 
by a contact (or point) transformation 

x = (f>(x, u, u x ), u = ip(x,u,u x ), (6) 

where 

d(p [dip di/j\ dip ( d(p d(f)\ 

7^«x + ^~ = 7^ T^U X ' 



OU x \OU OX J OU x \OU OX J 

All these equations can be linearized by simple differential substitutions (see [8]). 

In Section 3 we show that all equations (1) from our list admit a parametrization of 
the form 

q x = K(q, q y , q yy ), q t = S(q, q y , q yy , q yyy ), (7) 

where q x — K is a linearizable equation of second order and qt = S is a higher symmetry 
of this equation. For example, the parametrization (7) of equation (2) is given by 

„ _ %y _ _ %yy , q %y 
Qx — —5-, qt — ^- + 6 ^r- 

Qy Q 3 y Qy 
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Any common solution q(x, t, y) of system (7) gives us an one-parameter family of 
solutions for the integrable equation (1). However, the general solution of (7) depends on 
one function qo(y) = q(0, 0, y) of one variable, whereas the general solution of (1) depends 
on two functions of one variable. 

Although a general idea of such parametrization is contained in [9, 1], explicit formulas 
for linearization of equation (2) were obtained first by V. E. Adler. 
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2 Classification of third order equations. 
2.1 Integrability conditions. 

For equation (4) all mixed derivatives of q containing at least two time differentiation can 
be expressed in terms of 

Q, Ixi Qxxi ■ ■ ■ ■ ■ ■ ; Qti Qn — Qtx-, Qt2 — Qtxx, ■ ■ ■ , Qtii ■ ■ ■ (8) 

in virtue of (4). The derivatives (8) are regarded as independent variables. 
An equation 

q T = G(q, qi,q 2 ,..., q r , Qt, Qa, Qt2, • • • , Qts) (9) 

compatible with (4) is called infinitesimal (local) symmetry of (4). Compatibility implies 
that the function G satisfies the equation F(G) = 0, where 

F = & t -±^Dl-(£^D{) D t ^Dl-{M + ND t ) 

i=0 0< li \i=0 a Qti ) 

is the linearization operator for equation (4). 

In order to rewrite consistency conditions of (4) and (9) in terms of a series of conser- 
vation laws 

(Pi) t = (<n)x (io) 

for (4), one can use a formal Lax representation of the problem. The linearization of 
equations (4), (9) gives rise to the compatibility problem for linear equations 

cf> tt = (M + ND t ) 0, <f> T = (A + BD t ) 

or, equivalently, 

f 4> \ „ f 1 \ 



$ t = F*<f>, $ r = G*<S>, $ 
where 



4> t r \mn 



G *=(^|)' A d ^A t + BM, B^= B t + BN + A (11) 
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The cross differentiation yields 

D t (G*) = [F*,G*]+D T (F*) (12) 

where F#, are matrix differential operators. The crucial step in the symmetry approach 
(see [3, 5, 10] and references there) is to consider instead of above equation one as follows 

D t (R) = [F„R], (13) 

where R is matrix pseudo-differential operator. We call R matrix formal recursion oper- 
ator. 

Denoting as before R n = A, R 12 = B we can rewrite (13) as follows 

A tt - NA t + [A, M] + (2B t + [B, N\)M + BM t = (14) 

B tt + 2A t + [B, M] + [A, N] + ([B, N] + 2B t )N + BN t - NB t = 0. (15) 

If R\, R 2 are formal matrix recursion operators, then i? 3 = R\R 2 is a formal recursion 
operator as well and we find using (11) that 

A 3 — A X A 2 + B X B 2 M + B X A 2 ^ B 3 = A X B 2 + B X A 2 + B&N + B X B 2 ±. (16) 

Identities (14), (15) mean that the scalar pseudo-differential operator 1Z = A + B D t is 
related to the linearization T of equation (4) by 

J 7 (A + B D t ) = (A + BD t )F, (17) 

where A = A + 2B t + [B,N]. The pseudo-differential operator 1Z = A + B D t) whose 
components 

n m 

A = Y,a t Dl, B = Y,hDl 

— oo — oo 

satisfy (14) and (15), is called scalar formal recursion operator lor equation (4). If A and 
B are differential operators (or ratios of differential operators), condition (17) implies the 
fact that operator 1Z maps symmetries of equation (4) to symmetries. However we are 
using the notion of formal recursion operator for a completely different aim. 

Let 1Z\ — A\ + Bi D t and 1Z 2 = A 2 + B 2 D t be two scalar formal recursion operators. 
Then the product 71% = TZiTZ 2 , in which D\ is replaced by (M + ND t ) is also a scalar 
formal recursion operator whose components are given by (16). 

An operator S = P + QD t is said to be implectic if 

F*S + SF=0, S = P + QD t . 

Here and in the sequel the superscript * denotes the adjoint operator. If S can be applied 
to symmetries, then it maps symmetries to cosymmetries. In the symmetry approach P 
and Q are supposed to be formal non-commutative series with respect to D x . 

The operator equations of components of the formal implectic operator S = P + QD t 
have the following form 

P tt + U*P t + 2Q t V + QV t = V*P -PV- (QU + U*Q) V - U t *P, (18) 



4 



Q tt +2P t +2Q t U+U*Q t = V*Q-QV-(QU+U*Q) U -(PU + U* : P)-(U*Q+QU t ). (19) 

The linearization P + QD t of a variational derivative of any conserved density for equa- 
tion (4) satisfies equations (18), (19) up to a "small" rest (see [11]). 

Let us consider equations of the form (1). It follows from formulas (18), (19) that 
equation (1) has no higher conservation laws. Moreover, it is easy to prove that the 
density of any conservation law, up to total derivatives, is of the form 

p = r 1 (q,q x )q t + r 2 (q,q x ). 

All statements presented below can be easily reformulated for general equations (4). 
Theorem 1. If equation (1) possesses an infinite sequence of higher symmetries of 
the form 

Qn = Gi(q, q u q 2 , . . . , q u , q t , Qti, Qt2, • • • , q tSi ) (20) 
then there exists a formal recursion operator of the form 

K = (a + a^D' 1 + ...) + (D- 1 + b_ 2 D- 2 + . . .) D t , (21) 

where a,i,bi are some functions of the variables (8). 

For scalar evolution equations one can use (see [3, 5, 10]) the residues of powers of the 
formal recursion operator to derive the canonical conservation laws (10). Unfortunately 
for equations (1) this technique does not work and we present a different way (cf.[12]) to 
get necessary integrability conditions (10). 

Let TZ be a formal recursion operator of the form (21). It is then possible to find an 
operator 

TZ- 1 = (a^D' 1 + a_ 2 D~ 2 . . .) + (D~ 2 + /3_ 3 ZT 3 + . . .) D t (22) 

such that Tin- 1 = TZ~ l TZ = 1. Recall that we eliminate D 2 in virtue of T = in the 
product of scalar recursion operators. The operator is uniquely defined. 

Theorem 2. If TZ is a formal recursion operator of the form (21) for equation (1), 
then there is a unique representation of the total derivative operators D x and D t of the 
form 

2 3 

D x = Y jPi K\ D t = Y j a i K i . (23) 

— oo — oo 

Functions pi and <7j are densities and fluxes of some (maybe trivial) conservation laws (10) 
for equation (1). 

The next formulas define five integrability conditions (10) for equations (1): 

Pi = Ml, 
2 

p2 = v 2 + -a u 

1 o 

P3 = 6(T 2 - Mi(7i + 9u - 3wi1> 2 - 

O 

Pa = 6a"3 — 9ui<7 2 + 3af + 27u ui — u\ + 81t>i — 9u\v 2 — 27v 2 . 

p 5 = 2<t 4 + 18o"icr 2 — 27(<7i)t — 3>a 2 Ui — 3f7 3 iti — 9o" 1 -u 1 t> 2 — (J\u\ + 27aiu . 

The conditions mean that pi are densities of local conservation laws for equation (1). In 
other words, for any pi there exists a corresponding function o~i depending on variables (8). 
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2.2 List of integrable equations. 

Lemma 1. If equation (1) satisfies conditions (10) with i = 1,2,3,4 then it is of the form 

qtt = q X xx + {Mqt + A 2 ) q tx + A 3 q 2 xx + (A 4 q 2 + A 5 q t + A 6 ) q xx + 

(24) 

+A 7 qf + A 8 qf + A 9 q 2 + A 10 q t + A n , 

where the functions Ai depend on q and q x only. 

It is easy to verify that the class of equations (24) is invariant with respect to point 
transformations q — > (p(q). Moreover, if all functions A^ do not depend on q, then the 
shifts of the form q — > q + \\x + \ 2 t, where Aj are arbitrary constants, are also allowed. 

Theorem 3. Up to the transformations described above, any nonlinear equation (1) 
satisfying integrability conditions with i = 1,2, ... ,7, coincides with the equation 

qtt = q X xx + (3g s + k) q xt + (q t - 3q 2 x - 2kq x + 6V) q xx - 2V' q t + 6V' q 2 x + (V" + kV) q x , (25) 

where V(q) is any solution of equation of the form 

V' 2 = 8V 3 + k 2 V 2 + c{P + c (26) 

or with the equation 

(Z<k 3 \ 1 2 (2q 2 3q t \ / 3 2 s . ^ 

qtt = q XX x+ I — + 2 ) Qxt q * x ~ \ ~ + 2~ I qxx + C2 \ qxqt + 2 Qx )' ( ' 

\_ J y Qx Qx J 

where X(q) = c 2 q + C\ and q are arbitrary constants. 

Remarks. Actually, any equation (27) can be reduced to the equation with X(q) = 
const or to the equation with X(q) = q. The integrability conditions 6 and 7 have quite 
long expressions and we do not present them in the paper. We have only used these 
conditions to prove that any equation (1) for which all conservation laws (10) are trivial 
is equivalent to a linear one. 

2.3 Recursion operators. 

In this Section we present a closed form of recursion operators for models (25) and (27). 
The existence of these recursion operators implies the fact that all integrability conditions 
with i > 1 are fulfilled for these equations. 

Let us consider equation (25). This equation has the only non-trivial conserved density 
given by 

p = q t -q 2 x + 2V{q). (28) 
The simplest higher symmetries of (25) are 

q T = qxx + 2q x q t - q 3 x - kq 2 x + AVq x 

and 

q T = qxt + q x qxx + q\ + (q 2 x + 2V) q t - q x - kq\ + 6Vq 2 x + (V + kV) q x 
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They are generated by the following recursion operator 

K = D + (q t - 2q 2 x - kq x + IV) + q x D~ x (D t + 2q xx + 2V') 

acting on the seed symmetries q x and q t . A direct calculation shows that this recursion 
operator satisfies (14) and (15). In the degenerate case V = 0, k = we have K = 71%, 
where 1Z is defined by (3). The operator D t + 2q xx + 2V' corresponds to the variational 
derivative of the function (28) and therefore if we apply 1Z to any local symmetry admitting 
the conservation law with density (28), the result should be local. 
Another recursion operator for (25) has the form 

S = D t + (q xx -ql- kq 2 x + QVq x + kV + V) + qtD~\D t + 2q xx + 2V') 

One can verify that 

s 2 = n 3 -kns --n-c . 

2 

For equation (27) a recursion operator is given by 

ft = I \--X\- q x D I —D t + — h 




9x ^ 2/ 
The non-trivial conserved density for this equation is given by 

p = ^+ 1 -X(q). (29) 
Qx £ 

In the case X(q) — 0, Ci = the nonlocal variable Q = D~ l (p) satisfies equation (2). 

3 Linearization procedure. 

Each of equations (25) and (27) has only one non-trivial conserved local density. The 
Burgers equation u t = u xx + 2uu x possesses the same property: D t {u) = D x {u x + u 2 ) is 
the only conservation law for this equation. The crucial step in the linearization of the 
Burgers equation is to introduce the potential W of this conservation law. By definition 
the variable W satisfies conditions W x — u, Wt = u x + u 2 . It is easy to verify that the 
function U = exp {W) satisfies the heat equation Ut = U xx . 

For equations (25) and (27) the potentials of the conservation laws also play a key 
role in linearization. However the procedure of linearization is not so straightforward. To 
illustrate this procedure, let us consider the simplest version X = of equation (27). The 
potential of the conservation law for this equation satisfies the conditions 



w x = ± w t = q — + %. 

Qx Qx Q x 



A simple computation shows that the equation admits the following non-local symmetry 
q y = W . Since 




Ql 
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we have 

- w - -— - Si 

lyy — vv y — ■> %yy ~ 3 

or 

a — 1 a — qyVV 

ix — j qt — — — • 

Qyy %y 

After the Legendre transformation 

y = U z , q = U-zU z 
the latter pair of compatible equations becomes 

U x U zz , Uf U zzz . 
For the more complicated equation 



q X xx H q x t Q xx 5- q X x + c q xt q X x) (30 ) 

Qx q x qi v q x ' 



corresponding to the case X = const 7^ the same linearization scheme works also. The 
potential is defined by 

W x =*, W t = q -^ + q l + c*, 

q x q x qi q x 

the non-local symmetry is given by q y = exp (—cW). It is not hard to check that 

2 q y _ 3 QyQyyy ~ ^%%y 
q x — —c — , qt — c 3 . 

Qyy q yy 

After a contact transformation 

y = X - exp {-z) (U z + U), q= X - exp (z) (U z - U) 
these equations become 

2 3 
Ux = j (u zz - u), u t = c - (U zzz + U zz -U z - U). 

The most non-trivial case is X = q. In this case the potential W is defined by 

& 2 qr x q£ 2q x 4 

but non-local symmetries of the form q y = F(q, W) do not exist. However there exists a 
new non-local conservation law with potential Z defined by 

Z x = q 2 - 2Wq x , Z t = -^q 3 - 2q x - 2Wq t . 



Using these two potentials, we find a non-local symmetry q y = 2q exp (— j — ^-). Ex- 
pressing q x and q t in terms of the ?/- derivative, we get 

<l 3 <ly = Q 4 Qy(Q 2 Qyyy ~ ^QQyQyy + 

After the contact transformation y = z + jj-, q = — |j§ we do not obtain linear equations 
but 

To linearize the latter system one can introduce the potential Y such that Y z = U, Y x = 
-^j2 and after that make a point transformation Y <-> z. 

Equation (25) can be linearized as follows. It is easy to verify that it has a non-local 
symmetry q y = A(q) exp (-W), where 

W x = q t -q 2 x + 2V, W t = q xx - q 3 x + q t q x + k(q t -q 2 x + V)+ 6Vq x + V' + w 

and 

B ' 2 = B A + %B 2 + 8wB + b , B = -^ + \. 
It can be checked that 

q x = ^- + 2B(q), (31) 

Qy 

ft = + 3* + 3^3(g) + + 25(g)) - - B(q) 2 ) - (32) 

Function V from equation (25) is given by 

1 1 k 2 
V = -L B ' + -B 2 -—, 
4 4 48' 

where the parameters of the elliptic functions V and B are related by 

k A - 166 o 
ci = —32—, = 

After change of variables z = q,u = y equations (31), (32) take the following linear form 

u x = u zz - 2B(z)u z , 

3 \ k / \ k 2 

u t = -u zzz + 3B(z) u zz + - - B(z) 2 ) u z + -(w^ - 2B(z)u z ) + — u z . 

We see that in all cases there exists a non-local symmetry q y = G depending on the 
potentials such that q x and ft can be expressed in terms of y-derivatives by formulas 
(7). Using this parametrization, one can construct particular solutions of equations (25) 
and (27). 
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A parametrization of such kind arises not only for linearizable equations but also 
for equations of KdV-type and associated linear spectral problems. For example, let us 
consider the spectral problem 

V xx = {\ 3 + Ul \ 2 + u 2 \ + u 3 )V (33) 

For a non-local symmetry of this linear equation one can take (see [9, 1]) 

i|r , — — \U \T> — _L_ \T/ 

yy a 2 ' y/E 

Then 



1,1 /, , ^ 

where a(x, y) satisfies the Harry-Dym equation 



d X CI Ojyyy . 



At least on the local level the general solution of the latter equation depends on three 
arbitrary functions ao(x) = a(x,0), ai(x) = a y (x.O), a 2 (x) = a yy (x,0) of x and therefore 
this parametrization provides a generic potential in (33). Multi-phase solutions of the 
Harry-Dym equation lead to special potentials of the spectral problem (33). 
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